Periodic Modulation Effect on Self- Trapping of Two weakly coupled Bose-Einstein 

Condensates 
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With phase space analysis approach, we investigate thoroughly the self-trapping phenomenon 
for two weakly coupled Bose-Einstein condensates (BEC) in a symmetric double-well potential. We 
identify two kinds of self-trapping by their different relative phase behavior. With applying a periodic 
modulation on the energy bias of the system we find the occurrence of the self-trapping can be 
controlled, saying, the transition parameters can be adjusted effectively by the periodic modulation. 
Analytic expressions for the dependence of the transition parameters on the modulation parameters 
are derived for high and low frequency modulations. For an intermediate frequency modulation, 
we find the resonance between the periodic modulation and nonlinear Rabi oscillation dramatically 
affects the tunnelling dynamics and demonstrate many novel phenomena. Finally, we study the 
effects of many-body quantum fluctuation on self-trapping and discuss the possible experimental 
realization of the model. 
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I. INTRODUCTION 

Double well system is a paradigm model used to 
demonstrate marvellous quantum tunnelling property, 
e.g., cold atoms in optical lattice and magnetic spins in 
external magnetic field 0, ■ The realization of dilute 
Bose degenerate gas in 1995 provides a new opportunity 
to revisit the old topics 0, 0]. It should be addressed 
here, the tunnelling of Bose-Einstein condensates (BECs) 
system differs from the traditional quantum tunnelling 
mentioned above in two essential aspects: system scale is 
macroscopic (around 10/im), and more importantly it is a 
many-body system where the interaction between atoms 
plays an important role. A natural question arises that 
how the interaction between condensed atoms affect the 
quantum dynamics in double well system. This problem 
attracts much theoretical attention in past few years p| 
and recent realization of the BECs in optical double well 
trap bring new research surge 0, 0] • 

Among many findings, self-trapping is most interesting 
one |sl I9L lid flT| . It says, the BECs atoms in a symmet- 
ric double-well potential may shows highly asymmetric 
distribution as if most atoms are trapped in one well, 
even under a repulsive interaction between the degener- 
ate atoms. Recently, this counter-intuition phenomenon 
has been successfully observed in experiment Q- 

In this paper, with phase space analysis approach we 
study this phenomenon thoroughly. We emphasize the 
effects of a periodic modulation that is applied on the 
energy bias of the system. We find the external periodic 
field with high or low frequency can effectively modulate 
the transition parameters of the self-trapping. Moreover, 
for an intermediate frequency modulation, we find the 
BECs alternate between self-trapping regime and non- 
sclf-trapping regime with increasing atomic interaction 
parameter. This fact is attributed to chaos in phase space 
induced by the resonance between the periodic modu- 



lation and nonlinear Rabi oscillation. Furthermore, we 
consider the effects of many-body quantum fluctuation on 
self-trapping and find the quantum fluctuation will dra- 
matically influence the self-trapping phenomenon. Es- 
pecially for an intermediate frequency manipulation, the 
alternate phenomenon does not appear since there are no 
chaos in full-quantum description. 

The paper is organized as follows. Self-trapping ef- 
fects in two weakly coupled BECs is introduced and an- 
alyzed with phase space approach in section II. Section 
III discusses periodic modulation of this phenomenon for 
three cases: high frequency, low frequency and strong 
resonating region, respectively. After studying the effect 
of quantum fluctuation in section IV, we give the capable 
observation and application of our findings in section V. 
The acknowledgment is in the final part. 



II. SELF-TRAPPING IN TWO WEAKLY 
COUPLED BECS 



For two weakly coupled BECs, the wave function evolv- 
ing with time can be described with the superposition of 
individual wave function in each well. The spacial vari- 
able and the time one of the individual wave function 
are assumed to be separable. Under these approxima- 
tions, we can obtain a two-mode Schrodinger equation, 
i.e., two- mode Gross-Pitaevskii equation (GPE): 



dt 



= H 



(1) 



2 



where a and b are the probability amplitudes of atoms in 
two wells respectively. The Hamiltonian is given by 



H 



\b\ 2 - 



(2) 

where v is the coupling constant between the two con- 
densates in double-well potential, 7 is the energy bias 
between the two wells, c is the nonlinear parameter de- 
scribing the interaction 

[Hill El El- The total proba- 
bility is conserved and set to be 1. Our discussion focuses 
on the symmetric double-well trapped BECs system, so 
7 = 0. 

With a = \a\ e l9a , b = \b\ e l6b , the Schrodinger equation 
can be casted into a classical Hamiltonian system by in- 
troducing the population difference s = \b\ 2 — \a\ 2 and 
the relative phase = Ob — 6 a , 



H 



7s 



v\Jl — s 2 cos 8 



(3) 



with s, 8 are the canonically conjugate variables of the 
classical Hamiltonian system. 

Self-trapping motion refers to the trajectories whose 
average population difference is not zero < s >^ 0. This 
can be well understood in phase space of the classical 
Hamiltonian system (J3J). Three kinds of cases will emerge 
with different parameters, shown in Fig.l. 

1) For the weak interaction, i.e., c/v < 1, in the phase 
space, there are two fixed points p\, P2 at s = 0, 6 = tt 
and s = 0, 8 = 0, respectively (Fig. 1(a)). The atoms 
distribution corresponding to them is in equilibrium. The 
trajectories around them correspond to atoms oscillation 
between the two wells, i.e., Josephson oscillation, and 
the average of population difference to time is zero, i.e., 
(s) = 0. No self-trapping effect can be found. 

2) For the interaction 2 > c/v > 1, two more fixed 
points appear on the line 8 = tt. Among them, pi, p^ are 
stable andj»4 is unstable (Fig. 1(b)). They are at s = —k\, 



ki, 0, respectively, where k\ = y 1 — (f J- Obviously for 
fixed points p\ , P3 atoms distributions are not in equilib- 
rium. The average of population difference about the 
trajectories around them is non-zero, i.e., (s) 7^ 0. It 
indicates atoms self-trapped in one well. Because the 
population difference s and the relative phase 6 oscillate 
around fixed points, we denote it as oscillating phase self- 
trapping. 

3) With the interaction becoming stronger, i.e., - > 
2, new trajectories emerge, for example the trajectory 
across point p c (Fig. 1(c)). For these trajectories, s 
changes with time in [—1,0] or in [0, 1], relative phase 8 
varies monotonously. Obviously (s) ^ 0, atoms are self- 
trapped in one well, called run ning phase self-trapping 
predicted in the reference 0, H Eq | and observed in ex- 
periment 0- 

The above changes on the topological structure of the 
phase space associates with the change of the profile of 
the physical energy. With considering the relative phase 
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FIG. 1: Trajectories on the phase space of the Hamiltonian 
system (3) (left column). In right column we plot the energy 
profiles for the relative phase 9 = (dashed) and 8 — n 
(solid), respectively 



0.0 



A -0.5 

(fl 

V 



-1.0 



2 4 6 8 10 12 14 

c/v 



FIG. 2: The average of population difference vies the param- 
eter c/v when the trajectory s = — 1,6 = evolves with time. 



being zero or tt, the energy depending on the population 
difference s can be derived from Eq.(3) and plotted in the 
right column of Fig.l for varied parameters. The transi- 
tion from case 1) to case 2) corresponds to the bifurcation 
of the energy profile of 8 — 7r: The energy function of 
single minimum bifurcates into the function of two min- 
imums; The low limit of the energy profile with 8 = is 
— c/2, and the energy of the saddle point P4 is —v, which 



3 



is the local maximum of the energy profile with 8 = it. 
The transition from case 2) to case 3) is characterized by 
the overlap of the energy regime of the two profiles. In 
this case the trajectory originated from s = — 1,0 = 
can not cross the energy hill peaked by the saddle point. 
It will be confined in the low phase plane, corresponding 
the "running phase" type self-trapping. 

From the above analysis on the energy function we can 
obtain a general criteria for the self-trapping trajectories, 
that is H(s,9) < —v. With it, the transition parameters 
of self-trapping for the trajectory with arbitrary initial 
value can be found. For example, in the situation of 
recent experiment 7] s = 0.5,8 = 0, the c/v = 15.0 is 
the critical point for the transition to the 'running phase' 
type self-trapping. For the state of initial condition s = 
— 1,8 = the transition parameter is c/v = 2.0 in our 
analysis, agrees with our numerical simulations as shown 
in Fig. 2. The scaling law near the transition point shows 
a logarithmic singularity. 



III. PERIODIC MODULATION OF 
SELF-TRAPPING 

Generally, we can control the behavior of the system by 
applying a periodic modulation. For example tunnelling 
probability can be controlled through adjusting the mod- 
ulation parameters for the linear case |l7l lisl ITii l20| . 
In this section, we will discuss how a periodic modu- 
lation affect the nonlinear self-trapping. Without los- 
ing generality, we assume that the modulation is applied 
on the energy bias with amplitude A and frequency u>, 
i.e., 7 = As'mujt. Then the nonlinear GPE still can be 
mapped into the time dependent Hamiltonian (3). We 
then can investigate the global property of the trajecto- 
ries with the phase space analysis method as in above sec- 
tion. Differently, in this case, the phase space of the time- 
dependent Hamiltonian (or Poincare section) is drawn by 
stroboscopic plotting of the trajectories at the moment 
of integer multiple times of the modulation period, i.e., 
2ir/ui. We will consider following three typical cases cor- 
responding to different modulation frequency regimes. 



A. High frequency modulation {u> ^> v) 



above which the self-trapping occurs changes dramati- 
cally with periodic modulation. Fig. 4(a), taking the tra- 
jectory s = — 1, 8 = as an example, shows the critical 
interaction is about 1.531 at u> = 100, A/u> — 1.0, v = 1.0. 
Fig. 4(b) shows phase diagram for the transition to the 
self-trapping. It says that the critical values C/V de- 
pend on the modulation parameters in the way of the 
Bessel function Jq(— )■ 



in 




FIG. 3: Phase space of the Hamiltonian system (3) with 7 = 
Asinuit at uj = 100.A/oj = 1.0, obtained by stroboscopic 
plotting the trajectories with period 2-n /uo. 
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When the system is added by a high frequency mod- 
ulation, the two kinds of self-trapping, i.e., oscillating or 
running phase, still exist as shown in Fig. 3. Here the 
phase space are stroboscopic plotting of the trajectories 
at the moment of integer multiple times of the modula- 
tion period. Compared with Fig.l, the fixed points are 
shifted to the leftside. This is because the phase space 
depends on the time moment when Poincare surface of 
section is made. For the different time moment, the 
topological structure of the Poincare section keeps con- 
sistent which guarantee the validity of our observation. 
Moreover, we find that the critical interaction parameter 



FIG. 4: (a) The average of population difference vies param- 
eter c/v when the trajectory s = — 1, 8 — evolves with time 
at u) = 100,^4/0; = 1.0. (b) Transition phase diagram for the 
occurrence of the self-trapping. The triangle describes the 
case at u — 100. The line describes Jo (A/ui) (lower line) and 
2 Jo (A/uj) (upper one). The situation at uj = 10, 1000 also 
agree with the theoretical prediction through our calculation. 

The above observation can be understood from follow- 
ing deduction. Using high frequency approximation, the 



4 



time-dependent system Q is equivalent to a stationary 
one. Let us make the transformation a = e l ^ COSUJt a , 
b = e~ l 2^ cosult b . Then we can obtain the following 
Schrodinger equation 



.da 
'~dt 



db 



dt 

Using formula 



% , (4) 



(5) 



±iz cos ujt 



n= — oo 



Jn (z) (±t) 



(6) 



where J n (z) is the Bessel function of nth order, and con- 
sidering the contribution of those higher order Bessel 
function to the integrals can be neglected 0, |2l| , the 
effective coupling constant becomes v = vJo(^). 

Immediately, we can obtain the critical values are 
| = Jo (tj) f° r oscillating phase self-trapping and ^ = 
2Jo (— ) for running phase self-trapping based on the re- 
sults in section II. This result is consistent with the nu- 
merical calculation. 



± (c 2 / 3 — v 2 / 3 ) 3 2 0l- When the modulation frequency 
is very small, the parameter 7 varies adiabatically with 
time 7 = Aw cos uit Co. If amplitude A does not ex- 
ceed the critical values, atoms should be self-trapped. 



2/3 , „2/3\ 3 / 2 



for self- 



This gives the condition c > (A 
trapping, which is well coincident to the numerical result 
shown in Fig. 5(b). Below the critical point, the collision 
between the two fixed points mean occurrence of the in- 
stability, which may leads to the irregular oscillation of 
the average population difference as seen in Fig. 5(a). 
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B. Low frequency modulation (ui <C v) 

For this case, it is not significant to identify the two 
kinds of self-trapping. We suppose initially all particles 
confined in one well with the relative phase of n, i.e., 
s = — 1,9 = it. Then for given parameters we check if 
the trajectory is self trapped, < s >^ 0. We find that the 
self-trapping still occurs for some regime of parameters, 
as shown in Fig. 5(a). There, we see, above a critical 
value of interaction which is a function of modulation 
amplitude and the energy gap, the average population 
difference jump from a irregular oscillation to a constant 
((s) = —1, self-trapping regime). 

To see what happened before and after the transition, 
we plot the phase space of the system with 7 chang- 
ing slowly, seeing Fig. 6. Below the critical point, i.e., 
c/v = 13.5, the fixed p± (where we started from) moves 
up smoothly, at a certain j c , it collides with another fixed 
point pi and gives birth to a new trajectory. Above the 
critical point, i.e., c/v = 13.6, the trajectory pi moves 
along the line 9 = it up and down smoothly, having no 
chance to collide with other fixed point. For the large c 
case, the shift of our trajectory away from the bottom 
line (s = —1) is small, so that the time average of the 
population difference is approximately — 1. 

The singularity point 7 C where two fixed points col- 
lide leading to the topological change on phase space 
depends on the modulation amplitude but little on 
its frequency. Our theoretical deduction gives 7 = 



FIG. 5: (a) The average of population difference vies the in- 
teraction parameter at ui = 0.01, j4 = 10, v = 1.0. (b) The 
relationship between the critical values of c and A when u is 
extremely small at v = 1.0. The triangle describes the situa- 
tion at uj = 0.01. The line describes the theoretical prediction, 

c= [A^+v^y 12 . 



C. Strong resonating region (lj ~ v) 

For the strong resonating region, the situation is com- 
plicated. Due to the strong resonance between the 
nonlinear Rabi oscillation and the periodic modulation, 
chaos appears as shown in the phase space plotting Fig. 7, 
where the scattered points denote the chaotic trajectories 
forming chaotic sea[23|]. Inside the chaotic sea there are 
stable islands corresponding to the self-trapping trajec- 
tory. We find with increasing the interaction parameter, 
the stable islands appear alternately. Fig. 8 describes this 
alternate with sign = 1 standing for self-trapping and 
sign = for non-self-trapping when interaction parame- 
ter varies continuously. The mutations are the transition 
parameters. Analytic expression for the dependence of 
the transition parameters on the modulation parameters 
can not be found. 
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FIG. 6: Evolution of the phase space motions of the Hamilto- 
nian system (3) with 7 = Asinujt at A — 10, c/v = 13.5 (left 
collum) and A = 10, c/v — 13.6 (right collum) as 7 changes 
adiabatically. The arrows refer to the moving direction of the 
fixed points as (f> increases, i.e., 7 changes adiabatically, where 

4> = u>t. 



IV. THE EFFECT OF QUANTUM 
FLUCTUATION 

In the above discussions on the self-trapping, our 
framework is the mean field treatment with assuming 
that the number of particle is large enough. However, 
in practical experiment, the particle number is finite, we 
want to know how the quantum fluctuation caused by 
finite number of atoms affects the self-trapping. With 
this motivation we investigate following second quantized 
Hamiltonian which is the quantum counterpart of our 
mean field system (1), 

H = ~ (a f a - b]b )-^ (aWaa + tftfbb) + ^ (a^b + bh 

(7) 

where 7,c,v are energy bias between the two modes, 
the interaction between atoms and the coupling be- 
tween the modes, respectively. a\ (a, b) are cre- 
ating (annihilating) boson operators. N is the to- 
tal number of atoms. Its matrix elements arc ob- 
tained in the representation of Fock states \n,N — n). 
H nin = (N -n,n\H\n,N -n) ,.ff n ,„_i = if n -i,n = 
(N — n, n\ H \n — 1, N — n + 1). Corresponding to the 
mean field s = — 1 our initial state in the full quan- 
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FIG. 7: Phase space of the Hamiltonian system (3) with 7 = 
AsinLOt at ui = 2.0, obtained by stroboscopic plotting the 
trajectories with period 2n/cu. (a) describes the situation at 
c = 5.88, (b) at c = 9.2, and (c) at c = 16.65 respectively. 
The small fig is the magnified portion in the proper big Fig.. 



0.0 



A -0.5 

(/) 

V 

-1.0 



(a) 



1.0 



10 

C 



c 

.50.5 
n 



0.0 



(b) 



20 5 10 15 20 



FIG. 8: (a) The average of population difference vies the in- 
teraction parameter at lo = 2.0,^4 = 10, v = 1.0. (b) Critical 
values of c for occurrence of self-trapping shows complicated 
behavior at v = 1.0, w = 2.0, A = 10. See the text for detailed 
explanation. 



turn framework is \0,N). We then can solve the cor- 
responding Schrodinger equation i^a„ = H n;n -\a n -\ + 
,n+i a n+i with Runge-Kutta method, and 
trace the time evolution of population difference s = 
^2 l a "l ^~ 2n \ where a n is the amplitude of Fock state 
\n, N — n >. Our main results are that, for the cases 
of high or low frequency modulation, the self-trapping 
effect can still be observed in full-quantum description, 
i.e., with increasing the interaction parameter, the quan- 
tum average population will transition from symmetric 
distribution to asymmetric distribution. However, dif- 
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ferent from the mean field case, the transition here is 
smooth or continuous, no scaling law or singularity is ob- 
served. For the case of intermediate modulation, because 
the quantum effect dramatically suppress the classical 
chaos, the alternate phenomenon observed in the mean 
field case vanishes. What we observe is that a continu- 
ous transition to the self-trapping regime with increasing 
the interaction parameter. In following we give detailed 
presentation of our results. 

For to 3> v, with increasing the interaction parame- 
ter, Fig. 9 compares the transition to self-trapping in full- 
quantum treatment with that in mean field framework, 
by tracing the time evolution of population difference. 
It is found that, within mean field framework (Fig. 9(a)), 
there is an abrupt transition to self-trapping character- 
ized by a logarithmic scaling law as mentioned above. 
However, for full-quantum calculations ( Fig. 9(b)) the 
transition become smooth and no scaling law or singu- 
larity is found in this case. 
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FIG. 10: (Color online.) Evolution of population difference 
with time at u> = 0.01, A = 10, v = 1.0. (b) is mean-field 
results, (b ) is full-quantum results. In (a) the amplitude of 
line c = 14.0 is so small that it almost can not be seen. 
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FIG. 9: (Color online.) Evolution of population difference 
with time at ui = 100, A/uj = 1.0, v = 1.0. (a) is mean-field 
results, (a ) is full-quantum results. 

For case of w < w, our results are similar. Seeing 
Fig. 10, we find the transition to self-trapping regime is 
abrupt in mean field treatment (Fig. 10(a)), but is smooth 
for the full-quantum case (Fig. 10(b)). 

For the intermediate case of lo w v, the situation is 
quite different. Classical chaos shown in the mean field 
treatment is expected to be suppressed by the quan- 
tum fluctuation, alternate phenomenon observed in the 
mean field case vanishes, as we see in Fig.ll. In mean 
field method (Fig. 11(a)), self-trapping effect, (s) ^ 0, 
and non-self-trapping effect, (s) — appear alternately 
with increasing the interaction parameter, while, in full- 
quantum method (Fig. 11(b)), a continuous transition to 
the self-trapping regime is observed, that is, with increas- 
ing the interaction parameter, the BEC distribution be- 
come more and more self-trapped in a well. 

From the above analysis, we see in practical experiment 
to observe the dramatic effects of the quantum fluctua- 



FIG. 11: (Color online.) Evolution of population difference 
with time at u) — 2.0, A — 10, v = 1.0. (c) is mean-field 
results, (c ) is full-quantum results. 

tion caused by the finite number of particle, we need 
choose the parameters in the regime of the occurrence of 
classical chaos in mean field systems. 

V. CONCLUSIONS AND DISCUSSIONS 

We have investigated the periodic modulation effect 
on the self-trapping of two weakly coupled Bose-Einstein 
Condensates in a symmetric double-well potential. Ana- 
lytic expressions for the dependence of the transition pa- 
rameters on the modulation parameters are derived for 
high and low frequency modulations. It is shown that 
the self-trapping phenomenon can be controlled, saying, 
the transition parameters can be adjusted effectively by 
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the periodic modulation. We also study the effects of 
many-body quantum fluctuation on self-trapping due to 
finite number of the particle. We find that the quantum 
fluctuation make the critical phenomenon hazy for high 
and low frequency manipulation. For an intermediate 
frequency, the alternative phenomenon does not appear 
since there are no chaos in full-quantum description. 

Experimentally, BEC in optical double well system can 
be generated by coherently splitting BEC into a double- 
well potential generated by two laser beams. The radial 
separation of the potential wells is controlled by the fre- 
quency difference. The energy offset can be periodically 
modulated by adjusting the intensity of the two laser 
beams0,0|- The second possible physical system can be 
used to realize our model is the tunnelling between the 
internal state of BEC like 87 Rb[2J. There, two inter- 
nal states are separated by the relatively large hyperfine 
energy, but in presence of a near-resonant coupling field 



the states appear to be nearly degenerate. The energy 
bias can be adjusted by the detuning of the lasers from 
resonance. Our theory predicts the periodic modulation 
on the energy bias will dramatically affect the tunnelling 
dynamics and the nonlinear self-trapping. We hope our 
discussion will stimulate the experiments in the direction. 
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